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Abstract 

In this note, by making use of a hypergeometric series identity derived by 
Guillcra, I prove a Ramanujan-type series for the Catalan's constant. The con- 
vergence rate of this central binomial series representation surpasses those of all 
known similar series, including a classical formula by Ramanujan and a recent 
formula by Lupas. Interestingly, this suggests that an Apery-like irrationality 
proof could be found for this constant. 
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1. Introduction 

Catalan's constant, so named in honor to Eugene C. Catalan (1814-1894), 
who first developed series and definite integrals representations for it, is a clas- 
sical mathematical constant which may be defined as Q 

00 <—Y\ n 

G:=V-^ — '— = 0.9159655941 .. . (1) 
^ (2n + l) 2 

This constant is a special value of some important functions such as the Dirich- 
let's beta function /3(s) := J2n=o / ( 2n + ^ > namely /3(2) = G, and the 

Clausen's function 0^(0) :— 3(Li2(e* 9 )), nameljQ 



CI: 



and Cl2(37r/2) = — G. For positive integer values of n, we can trace an anal- 
ogy between j3{n) and C,(n) := 2fcLi V^™) ^ ne Riemann zeta function, since 
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of z, z < 1, and extended to C, except for the cut [1, oo), by analytic continuation. 
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both C(2n) and (3(2n — 1) are rational multiples of 7r 2 ™ and 7r 2 ™ -1 , respectively, 
whereas finite closed-form expressions for both ((2n + 1) and /3(2n) in terms 
of other basic constants are unknown However, the proof by Apery (1978) 
that C(3) is irrational Q has created an 'asymmetry' in that analogy because 
the irrationality of /3(2), though very suspected, remains unproveno 

From the point of view of numerical computation, Catalan himself (1865) 
computed G to 14 decimal places 0]. By making use of a technique from Rum- 
mer, Bresse (1867) computed it to 24 decimals, a result that was improved to 
32 decimals by Glaisher (1913) Q. With the advent of computers, G has been 
computed to a large number of digits. For instance, Yee and Chan (2009) com- 
puted it to 31 billion decimals [15j. Their computation employs two formulas, 
one of which is a central binomial formula due to Ramanujan (1915) (To]|FI 

G = - In (2 + V3) + - V -r- . (3) 

On searching for similar rapidly converging series, Lupas (2000) has found the 
following alternating series [9j 

64 h n3(2n-l)( 2 :)(-) 2 

This series converges so fast that it has been implemented in Mathematica™ 
(version 6) for computing G. 

On searching for new congruences modulo primes, Z.-W. Sun (2011) has 
pointed out that the following central binomial series should converge to G (see 
Conjecture A7 of Ref. [H): 

He indeed comments that this formula could be derived from a hypergeometric 
identity proved by Guillera in a recent work [6|, but a complete proof is not 



provided there in Ref. [13 1 



Here in this note, starting from a Guillera's hypergeometric identity I prove 
a Ramanujan-type series representation for Catalan's constant similar to that 
in Eq. ([5]) whose convergence rate surpasses that of all known central binomial 
series representations for this constant. 



2 Presently, the only known irrationality results for even beta values are the recent proofs by 
Rivoal and Zudilin (2003) that there exist infinitely many positive integers n for which /3(2n) 
is irrational, and that at least one of the seven numbers /3(2), . . . ,/3(14) is irrational [Till . 

3 This can be proved from the fact that G = — Jq^ 4 In (tan 9) dd = — | f^ 12 In (tan 9) d9, 
as nicely described in Ref. Q. 
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2. A new Ramanujan-type series for G 

Let us adopt the usual notation for the generalized hypergeometric series: 
01, . . . , a p \ _^ (ai) n ■ • • (%>)„ z 



pF AC'^b/' z )-^ o ( b S n ...(b q ) n n nr (6) 

where (a) = T(a + n)/r(a) is the Pochhammer symbol. Our main result 
makes use of the lemma below, which determines a special value for 3F2 f ai,a ^ ,a j i 

\ 01,02 

a function that converges at z = 1 whenever 5i{(oi + 62) — (a% + a 2 + 03)} > 
(see, e.g., Eq. (2.2.1) of Ref. 0). 



Lemma 1 (A special value). 

3F 2 ( H'l \\]=2G. 

V 2' 2 / 

Proof. We start from a well-known integral representation for generalized hy- 
pergeometric functions (see, e.g., Eq. (1.2) of Ref [7|), namely 



(7) 

valid whenever 3?(a) > and 9? (7 — a) > 0. It then follows that 



2' 2 / M^M'J « v 2 



1 ,.1 2i*l I 3 i 



2 Jo 



2*1 1 h j ;x 2 ) dx. (8) 



'0 V 2 

Now, let us show that, for all x g (0, 1), 



2 



. 1 1 2 \ arcsinz 
2 FW ' 3 ;x 2 = 7 ==- (9) 



It is well-known that 

2*1 (H IX 2 ) (10) 

for all non-null values of x for which th e hy pergeometric series at the left-hand 



side converges (see Eq. (1.5.10) of Ref. [12|). Two sucessive applications of the 
Euler transformation formula 



3 



on Eq. (fTU)) lead us to 



arcsm x 

x 

and, after some algebra, 

arcsin x 



1 - x 2 



2*1 



3 i 



1-X* 



2 F 1 \ H -x 2 



X vi - X 2 \ 2 

This completes the proof of Eq. ©. From Eqs. ((HJ and ©, one has 



5=1,1 _ i 

/ JO •'• vT ,r- 



3^2 ( 2 ' 3 ^ ;i 



The trigonometric substitution x = sin 9 reduces this integral to 



1,1,1 

3 3, 
2' 2 



3*2 1 2 '.^ ;1 I = 



7T/2 g 



sine) 



d<9, 



(12) 



(13) 



(14) 



(15) 



which can be solved in terms of the dilogarithm function Li 2 (z) as follows. 
First, note that 



dO 



In 



1 — exp {it 



+ i[U 2 (-e ie )-U 2 (e w )], (16) 



1 + exp(i6>)_ 

as can be easily checked by differentiating the right-hand side. Then 
•*/a e 



dO 



o 



suit 



7T / 1 — 2 

— In ; I + i [Li 2 (— i) — Li 2 (i)l — { lim a In 

2 V 1 + 1 J 1 a^0+ 



1-el 



1 + e( iQ ) 



+ i[Li 2 (-l)-Li 2 (l)] 



7T 



hi 



-2<Cls 



+ i 



7T 7T 

12 ~ "6 



. | M _ 0+ja .(:) +4 i.-f ( ta i-.f) +M+ «V 

= 2G, 

where the special value of the Clausen function in Eq. ([2]) and the principal value 
of the logarithm function, with Arg(z) £ (— 7r, tt], were taken into accountQ The 
substitution of this result in Eq. (fTS"]) completes the proof. 



(17) 



□ 



4 On Entry 9 of Adamchik's webpage lj, where several representations for G are proved 
computationally with Mathematica™ , one finds the integral formula i Jq*^ 2 Q I sm S d6 = G. 
Our Eqs. OH and (T?) can then be viewed as a formal proof of this formula. 
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We are now in a position to prove a rapidly converging central binomial 
formula for the Catalan's constant, which is our main result. 



Theorem 1 (Rapidly converging central binomial series for G). 



C - 1 yv 1)n (3n + 2)8» 
2^ l (2n+l)3( 2 ;) 



Proof. Let 



/(*) := £ (-1)" g iy + ^3 [6(* +») + !] (18) 



n=0 

be a function of a real variable a; in the open interval (0, 1). Then 

) 3 (6n + 4) ^ n! 3 (6n + 4) 



= 2 n! 6 (3n + 2)64" 



OC 



On the other hand, from Guillera's third identity in Ref. @, we know that 
/ W -4«f ^ + ».^ + ». (20, 
for all values of a; for which this series converges. Therefore 

/ 2x + l 2x+3 i \ 

m=^ 3 F^ (21) 

which implies that 

/(T) =2 3 F 2 fK 1, 3 ;lV (22) 



2' 2 



ill 

On substituting the result obtained in Lemma Q] for 3F2 ( 2 '3 '3 ; 1 ) , we find 



/(|) = AG. From Eq. (fl9j). one then has 

^ (2n + l)3(^) 3 

which completes the proof. 



2' 2 



(23) 



□ 



The convergence rate of the just proved central binomial series representation 
for G is to be compared to that of other known similar series, including those 
in Eqs. ((3|) and (j4]). This is done in details in the next section. 
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3. Convergence rates 

Let us now check the convergence rate of each Ramanujan-type series for G 
mentioned in this work. By applying the Stirling's improved formula, namely 



n!~(3f-( n+ £)' (24) 

which implies that 



2n\ 2 2 « v /2n+i 



we shall develop an order estimate of the n-th term for each series. 



(25) 



We begin with Ramanujan's series, given in Eq. ([3]). Its n-th term is 
1 



(2n+l) 2 ( 2 n n ) 2 2n + ± (2n + l) 2 

6n + 1 6n + 1 



3 2 2n y/12n + 1 (2n + l) 2 2 2n V^n + 1 (2n + l) 2 
3 



2 2 ™ (2n + 1) y/12n + 1 



(26) 



The more complex central binomial series by Lupas, see Eq. has an n-th 
term whose absolute value can be estimated as follows: 

2 8n |40n 2 - 24n + 3| 2 8 ™ |40n 2 - 24n + 3| 



n 3 (2n - 1) ( 2n ) (i n ) n 3 (2n p 22 " V 2 "^ 2 8 " 3 24 " +1 .. 

|40n 2 - 24n + 3| (12n + l) 2 (n + ±) 



n 3 (2n - 1) 2 2n ^2n + | (24n + 1) 
For sufficiently large values of n, this simplifies to 
(40n - 24) (12n + l) 2 v 7 ^ (n + ±) (40n - 24) (12n + l) 2 (6n + 1) 



(27) 



2 2„ n 2 (2 „ _ 1} ^2^71 (24n + 1) 2 2 " ™ 2 ( 2 ™ - 1) V^ + T (24n + 1) 

4 (5n - 3) (12n + l) 2 (6n + 1) 4 (5n - 3) (12n + 1) (6n + 1) § 

" 2 2 " n 2 (2n - 1) Vl2n + 1 (12n + ±) 2 2 ™ n 2 (5n - §) Vl2n + 1 _ 

10 (12n + 1) (6n + 1) _ 10 v / 12n+ 1 (6n+ 1) 

2 2 " n 2 v / T2nTT ~~ 2 2n n 2 

_ 5Vl2ir+T(6n+l) 

• ( 28 ) 
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This convergence rate is clearly slower than that of Ramanujan's series, not to 
say the number of basic arithmetic operations needed to compute the n-th term, 
which is considerably larger in the Lupas' seriesH 



The convergence rate of the series we proved in our Theorem [T] can be esti- 
mated as follows. The n-th term of our series is 

(3n + 2) 2 3n (3n + 2) 2 3 ™ 



(2n+l) 3 ( 2 „") d (2n 



2 6n+3 



(12n+2)7 

_ 3(n+ §)2 3 "(12n + 2)i 
8 (n + ±) 3 26«+3 

For sufficiently large values of n, this can be approximated by 

3 (12n + 2)i _ 3^/2 (6n+l)* 



(29) 



8 (n + i)2 23«+3 (2n + l) 2 2 3n + 3 

9 V2 7 ^S^^27V2- 



(2n + l)23«+3 (V6^+ 3) 2 2 3 «+ 3 

27V2 - 1 5/V * (30) 

V6n + 3 2 3n + 3 2 3 " y/2n + 1 

The factor 2 3n makes our series converges faster than both the Ramanujan and 
Lupas' series. 

In Table [l] below, we compare the error committed in approximating G by 
the partial sum of the first N terms corresponding to each central binomial 
series mentioned in this work. It is clear from this table that our series is 
the one that yields the smallest absolute error. The only competitive series is 
that conjectured by Sun, see our Eq. ([5]), but a direct comparison of its n-th 
term with that of the series in our Theorem [T] shows that it converges slower. 
Therefore, even an eventual proof of Sun's conjecture will not furnish a central 
binomial series faster than the one I am presenting here in this work. 
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5 Despite these disadvantages, Lupas' series has been implemented in Mathematica™ (ver- 
sion 6) for computing G [T3jl . 
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Tables 



N 


Lupas 


Ramanujan 


Sun 


Theorem [T] 


5 


+2.9 x 10" 


-4 


-3.0 x 10- 


-6 


+1.1 x 10" 


-5 


-1.3 x 10~ 6 


10 


-2.0 x 10" 


-7 


-1.3 x 10- 


-9 


-2.6 x 10" 


-10 


+3.1 x 10~ n 


50 


-7.7 x 10- 


-32 


-1.1 x 10- 


-34 


-9.1 x 10" 


-47 


+1.1 x 10" 47 


100 


-4.3 x 10" 


62 


-3.3 x 10- 


-65 


-4.5 x 10" 


-92 


+5.6 x 10~ 93 


500 


-2.9 x 10" 


-303 


-4.6 x 10" 


-307 


-1.2 x 10" 


-453 


+1.5 x 10~ 454 


1000 


-1.9 x 10" 


-604 


-1.5 x 10" 


-608 


-2.4 x 10" 


-905 


+2.9 x 10~ 906 



Table 1: Deviations from G of the partial sums obtained by adding the first N terms of each 
central binomial series mentioned in the text. 
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